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I. INTRODUCTION 



In perturbative double Compton scattering [HI2JE!, an incoming photon interacts with an 
electron, and two photons are emitted. This process, which is represented by the Feynman 
diagrams in Fig. [TJ can be described by perturbative quantum electrodynamics (QED) and 
requires no other special theoretical input. Experimental evidence ranges from the first 
measurements more than 50 years ago [U [6] to the more recent [TJ [9J [10] • However, if 
the emission process takes place inside an intense laser field, then the physics changes, and 
the electron line is dressed by multiple interactions with the laser field (see Fig. [2j. The 
emission of two photons is a purely quantum process which cannot be described by classical 
radiation theory [TT]. An exception is encountered only for the case of the sequential emission 
of two quanta which occurs when the intermediate propagator hits a resonance pole, given 
by a resonant Dirac-Volkov state. In that case, to which we will return to later in the paper, 
the diagrams in Fig. [2] break apart into two distinctive blocks for the emission of the two 
photons. 




FIG. 1: The 6 Feynman diagrams contributing to (ordinary) double Compton scattering (initial 
electron four momentum pi, final four momentum pf). The electron line is denoted by a single 
customary fermion line. The frequencies of the two emitted photons with wave four vectors kb and 
k c may be different and thus may be their wavelengths; this is explicitly indicated in the panel. 

The most interesting geometry for the process is the backscattering case, where a relativis- 
tic electron counterpropagates against an intense laser beam of comparatively low frequency 
(on the order of a few eV) . In ordinary Compton scattering, the electron is usually assumed 
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FIG. 2: Feynman diagram for the two-photon decay of a Dirac-Volkov state. The electron line is 
dressed by the laser field and denoted by a zigzag laser photon superimposed on the fermion line. 

to be at rest, and the scattering of a highly energetic photon is considered. Because the 
kinematics is inverted in the backscattering case, one sometimes refers to this scenario as 
"inverse Compton scattering" . During the emission, the electron interacts with the laser field 
via an arbitrary number of interactions (see Fig. [2J; the process can be described by fully 
laser-dressed Dirac-Volkov propagators p21 HB]- So, we may refer to the process depicted 
in Fig. [2] as "inverse laser-dressed double Compton backscattering." 

Note that for a single-photon Compton backscattering, the highest photon energy at- 
tainable is A^uj, where us is the laser photon energy and 7 is the Lorentz factor of the 
incoming electron. For a defined scattering geometry, the energy of the emitted photon 
thus is uniquely defined, and it coincides with the energy of the emitted classical (Larmor) 
radiation in the specified direction provided the laser photon energy is much smaller than 
the electron mass and Lorentz boost factors are taken into account. If the electron absorbs n 
laser photons during laser-dressed single Compton scattering, the energy maximum changes 
to An^ 2 uj/{1 + £ 2 ), where the laser intensity parameter £ is defined in Eq. ^ below (£ 2 is 
proportional to the laser intensity). When two photons are emitted in laser-dressed double 
Compton scattering, their maximum energy sum is limited by U\ + uo 2 < 4n7 2 u;/(l + £ 2 ). 
As we will show, it is possible to designate energy and angular regions in which the dou- 
ble scattering process dominates over single scattering, which is crucial for an experimental 
verification [HI [15] . 

Interestingly, as noted in [161 EL71 IS]) the two photons emitted during the process are 
entangled because of the quantum nature of the process. In order to quantify the en- 
tanglement, the emission directions of the two quanta cannot be used with good effect, 
because they represent continuous variables in three dimensions. However, the polarization 
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components of the two photons along the emission lines can be uniquely decomposed in a 
two-dimensional space composed of unit vectors (effectively a one-dimensional space), and 
measured independently Triggering on simultaneous two-photon events, one can then mea- 
sure the entanglement quantitatively: an appropriate measure is the so-called concurrence 
[19] 120] which measures the polarization entanglement of the two quanta. 

The usual double Compton scattering, which involves the absorption of only one laser 
photon, has a rate which is proportional to the square of the laser four- vector amplitude, i.e., 
proportional to its intensity. Therefore, we may refer to the single scattering process as the 
"linear" process. With rising laser intensity, the rate deviates from the simple linear intensity 
dependence, it becomes more and more indispensable to include higher-order effects, and 
the process becomes nonlinear. 

In order to bring the current investigation into perspective, we would like to mention other 
work performed in connection with two-photon emission from free electrons: indeed, a pair 
of photons may be produced by electrons accelerated by any kind of external field. Probably 
the most well known process of this kind is double bremsstrahlung [2U [221 1231 [2U [23 1221 [23 
281 [29] , but also photon pair creation in a magnetic field [301 EH E2] , and in a crossed field 
[33] has been considered. The process under investigation in this paper is complementary to 
those mentioned above, and may provide for better control of the properties of the produced 
photons by adjusting the laser parameters. 

This paper is organized as follows. In Sec. [TT| we discuss the formulation in terms of 



a laser-dressed ( "nonperturbative" ) QED formalism. We then continue, in Sec. Ill, with a 
comparison of the predictions of the fully relativistic, nonperturbative theory to the relativis- 
tic, but perturbative (in the laser field) theory of double Compton scattering. In particular, 
we extend the discussion given in Ref. [31] to also include circularly polarized laser fields. 



In Sec. IV, we study the angular correlation and the entanglement of the emitted photons, 
in the nonperturbative formalism. Finally, conclusions are drawn in Sec. |Vj Throughout 
the paper, we use relativistic natural units such that h = c = 1, and a space-time metric 
gtiu _ diag(l, — 1^ — 1 ; —1). Scalar products of four-vectors are written as p^q^ = p- q for two 
four- vectors p and q. The gamma matrices are written as 7 M , and their contraction with a 
four-vector p as p = 7 • p. 
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II. FORMULATION OF THE QED THEORY 



A. Notation 

The electron mass is denoted by m, and the electron charge by e = — |e|. The laser wave 
vector points in the negative x 3 direction [with the space-time coordinate x^ = x) = 

(x°, x 1 , x 2 , x 3 )}, 

x*=(u,x)=u (1,0,0,-1), (1) 

and the laser four-vector potential, modeled as a monochromatic plane wave, for linear 
polarization is 

Ag n (0) = o"cos^ <P = x-x, (2) 
with a ■ k = 0, = \a\(0, 1, 0, 0). For circular polarization we have instead, 

^circ(^) = a l COS 4> + a 2 4>i (3) 

with ai • a 2 = 0, a\ = a\, = |«i | (0, 1, 0, 0), = |ai|(0, 0, 1, 0). The laser intensity 
parameter £ is defined as 



— e I —a 2 



£=— \— > ( 4 ) 



for linear, and 



m m V 2 



,2 „2 



for circular polarization. For a consistent comparison of linear and circular polarization, one 
should compare at the same value of £, which corresponds to the same laser intensity. The 
parameter £ relates to the root-mean-square electric field amplitude E like 

( = (6) 
muj 

and can be said to be the relativistic (inverse of the) Keldysh parameter: £ < 1 corresponds 
to the multiphoton regime of relativistic laser-matter interaction, where the coupling to 
the laser field is perturbative, and £ > 1 is commonly referred to as the tunneling, or 
nonperturbative regime. The quantum parameter x [35], which in general determines the 
magnitude of quantum effects such as e + e~ pair creation, spin effects etc, is defined as 

X = i — a - , (7) 
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where pi is the initial momentum of the electron [see Eq. (10)]. If we compute x i n the rest 
frame of the electron, where Pi = (m, 0), then 

u E 



(8) 



m E crit 

where E cv - lt = m 2 /\e\ is the critical (Schwinger) field. Thus, x is the amplitude of the 
electrical field of the laser compared to the critical field in the rest frame of the electron. 
The relation to laser intensities follows from the formula 



£ 2 ( ) ^crit; 



(9) 



\mJ 

where I CI a = 2.3 x 10 29 W/cm 2 is the critical intensity, corresponding to x = 1 m the lab 
frame. For most of our examples, we use to = 2.5 eV (optical laser), which corresponds to 
/ = 5.5 x 10 18 W/cm 2 for £ = 1 and / = 2.2 x 10 19 W/cm 2 for £ = 2; the laser field here is 
strong but manifestly sub-critical. Note that even in the case of a relativistic (Lorentz factor 
7 = 10 3 ) electron beam as considered later in the examples (see Sees. HE and III), the laser 
field remains sub-critical in the rest frame of the electron, since % = t^pi ■ x/m 2 w 10~ 2 1 
with the parameters chosen. 

The initial electron four-momentum is (we assume the electron to be counterpropagating 
with respect to the laser field, i.e., moving in the positive x 3 -direction) : 



Pi 



{E hPi ) = {E i M \fE\ 

,2 



m 



x 



(10) 



2k- pi 

which is valid for both circular and linear polarization. The final electron four-momentum 
is 

-2 



Pf= ( E f,Pf), Qf = Pf + C 



m 



2x • pf 



x= (Qf,Qf)- 



(11) 



The four- vector q it f introduced in Eqs. ( 10 ), ( 11 ) is the average momentum of a laser-dressed 
electron [35], with corresponding average mass m*, 



Qi 



m 2 (l + e 



(12) 



The electron spinors are used in the following form: 

/ 



U r (p) 



E + m 
2m 



5 rl 

5r2 



\ E+r 



■Pit)) 



(13) 



with the standard vector cr being composed of the (Pauli) 2x2 spin matrices. With this 
convention, the spinors are normalized according to u\.{p)^°u r (jp) = u r u r = 1. For an 
electron moving in the x 3 -direction, r = 1 corresponds to a right-handed electron, and r = 2 
to a left-handed electron. 

The Volkov states (35], solutions of the Dirac equation with an external laser field 



(id -m- eA^j * = 0, 
read for linear polarization [see Eq. ([2|] 




A {s,a,(3) + 6Xa Ai(s,a,/3) 
Ik ■ p 



(14) 



x u r (p)e-^ q+s ^ x , 



where 



a 



ea ■ p 



P 



e 2 a 2 



x ■ p 8x ■ p 

Here, the generalized Bessel function [351 EZ] is defined as 

A k (n,a,(3) = — 



(15) 



(16) 



cos kQ e ine-i aS ine+il3sin2e - 



(17) 



with k > 0, from which follows v4fc >0 (n, a, (3) = \ [Ak-i(n + 1, a, (3) + Ak-i(n — 1, a, (3)]. 
For circular polarization [see Eq. pj] we have, 




E 



. exa 2 

+ 2^~p JAa ^\ 



u r (p)e 



2x ■ p 



-i(q—sx)-x 



Here 



and 



2,2 eoi • p ea 2 • p 

a = \/af + a 2 , «i = , a 2 



x • p 



x • p 



<£> = arctan 2 (a 2 , —on) ■ 
The arctan 2 (■, •) functions is defined as 

arctan 2 (y,x) = arctan ^— j if x > 0, 
arctan 2 (y, x) = ix + arctan ^— j if x < 0, 



(18) 

(19) 
(20) 



(21) 



the usual Bessel functions are denoted by J n (a), and 



(22) 



Note the normalization factor in Eqs. (15) and (18): the volume V comes with the wave 
function, and not with the spinor u(p). 



e 




if, = tt/2 ip = 3tt/2 





FIG. 3: Polarization vectors according to Eq. (24) for fixed, small 6 ~ 0, i.e., for a photon 



propagating in the positive z direction. For example, we have at ip = ir/2 the two vectors 
e 1 = (0, 0, cos 0,- sin9) (0,0,1,0) and e 2 = (0,-1,0,0). 



The propagation four-vectors of the two emitted photons are denoted by 

k b = (uj b) k b ) = uj b k b 

= 0^(1, sin 6 b cos ip b , sin 9 b sin ip b , cos 9 b ), 

k c = (^C; kc) = &c 

= u c (l, sm9 c cosip c , sm9 c smip c ,cos9 c ), 



(23a) 



(23b) 



ip measuring the azimuth and 9 measuring the polar angle. As a basis for the two polarization 
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four-vectors e b and e c of the two emitted photons, we take 

4 = (0, cos6*;, cos ipb, cos 6^ sin-^b, — sin^b) , 

el = (0, -sin^ b , cos ip b , 0) 

el = (0, cos# c cos ip c , cos 9 C sin ip c , — sin# c ) , 

e\ = (0, -sinip c , cos^ c , 0) . (24) 

As an aid to the discussion, Fig. [3] illustrates the direction of the polarization vectors for 
small polar angle 9 and different values of ip. Alternatively, the polarization can be expressed 
in a helicity basis according to 




FIG. 4: Clarification of the index s. Shown above is one of the contributing Feynman diagrams in 
the perturbative picture, where the laser photons are inserted one by one. The net number of laser 
mode absorbed photons in this case is n = 1. The propagator momentum is p = qi — k b + sx, so 
that s counts the net number of absorbed photons before emitting photon k c , i.e. the momentum 
at the position of the label u p" . For the above diagram, s = 0. Although n must be positive for 
a net two-photon emission process, s may be negative, and to get the total amplitude for fixed n, 
one should sum all diagrams of this kind with s ranging from — oo to +oo. 

B. Matrix element for linear laser polarization 

The S'-matrix element for two-photon emission from a Dirac-Volkov state follows from 
standard Feynman rules, with four-vector potentials A b = A^-y^ = e b e lkb ' x / \JV 2u b and 
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A c = A^^ = e c e tkc ' x I \JV 2u c for the two emitted photons (see also the Feynman diagram 
in Fig. |2l). For linear laser polarization, we get 



olinear 



c(l) , c(2) 



(26) 



d Xxj d x 2 4> qjrf (x 2 ) iA c (x 2 )iG(x 2 ,x 1 )iA b (x 1 ) 

+ iA b (x 2 )iG(x 2 , Xi)iA c (xi) ip qun (xi) 
(2vr) 4 e 2 m 



^ 2V 2 y /u c u b Q i Q f 



8 (li - qj + nx - k b - k c ) 



x u rf (p f ) 



+M s r 



M s b ~ n 



2 m z 



2x-p b 



x + m 



2 



Pb ~ m l 

■x + m 



P 2 C - ml 



U r , {Pi 



Here, G(x,y) denotes the laser-dressed propagator function [131 EH], which can be con- 



structed from the Volkov state (15). The propagator momenta are given as 

q { + sk - kb, p c = qi + sk - k c . 



Pb 



(27) 



The matrix element is proportional to V' 2 , since there are one in-state and three out-states, 
each with a factor y/V. Here, n is the net number of absorbed laser photons, and the 
summation index s can be understood as the number of laser photons absorbed up to and 
immediately before emitting the second photon (see Fig. [4] for a pictorial explanation). 

The matrix-valued functions for the transition currents M and F are given as follows. 
For the first channel, we have 



M{;- n 



A (s -n,a f - a b ,/3f - P b ) e c 



+ A 1 (s-n,af-a b ,f3f-l3 b ) e 



2x ■ pb ^ 2x • pj 



6^0^% m 6 

- A 2 (s -n,a f - a bl p f - (3 b )- — x, 

2x • pfX ■ p b 



(28) 



and 



F b s = A (s, on - a b , pi - p b ) h 
+ At(s, a>i - a b , Pi - p b ) I e b 



+ 



2x • pi 2x • pb 



- A 2 (s, ai - a b , p { - p b ) x . 

2x ■ p b x ■ pi 



(29) 
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For the second channel, the two currents are given as M s c ~ n = M£~ n (c <-> b) and = 
F b s (c <-> b) under replacements of the corresponding expressions for the first channel. The 
arguments entering the generalized Bessel functions read 

ea-Vj _ e 2 a 2 . . 

with j G {i, /, 6, c}. The spinors it r . describe the spin state of the in- and outgoing electron, 
respectively. Note that x-Pij = X'Qi,f, and that due to x 2 = 0, a Cj b and /3 C)b are independent 
of the summation index s, although one might have initially assumed a dependence on s in 
view of the presence of p b and p c in their respective defining equations. 



C. Matrix element for circular laser polarization 



For the case of circular polarization of the laser, the matrix element can be derived in a 



similar way to Eq. (26). The matrix element reads 



ocircular c(l) c i q(2)c 

°fi — °fi ' °fi 



(31) 



(27r) 4 e 2 m 



2V 2 ^/uj c uj b Q i Q f 



$ 4 (<li - qf + nx - k b - k c 



x u rf (p f ) 



Pb - £ 2 t s —*c + m 

1V b „ 2 __2 G 6 



+N s c ~ n , * G s 

pi - m i 



U n (Pi 



Here, as is typical for circular polarization, the generalized Bessel functions in the formulas 
simplify to ordinary Bessel functions. The matrix-valued functions for the first channel read 



e 2 a 2 x ■ e c 



-K 



+ 



2x • pfX ■ pb 



(32) 



1 / A exai eaix 
x - I e r - h e, 



2 \ 2x • p;, 2x • pf 



1 / exa 2 ea2X 

x — e c 1 e< 

2z \ 2x • ^ 2x • 
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and 



1 2k ■ ■ p b ' 



1 / ekdi eh\k 

x o e b - + e b 

2 \ 2k ■ Pi 2k ■ p b 



1 / eka 2 ea 2 k A 
2% \ 2>c ■ pi 2k ■ p b 

For the second channel, we have N^~ n = N b s ~ n (b <-> c) and G s c = G s b (b <-» c). Here, 

5 /6 = v Kfc) 2 + («y 2 > «/c = v Kc) 2 + Kc) 2 , 

x _ epf ■ a x ep b ■ a x 2 _ e Pf ' a 2 e Pb • a 2 

a fb ~ > a fb — ) 

J X'Pf X ■ Pb K 'Pf H '■ ' Pb 



(34) 



1 2 

and similarly for a^c, a ib \. The phases <£> can be expressed in terms of the generalized 



arctan function (21) 



<p fb;C = arctan 2 (-a 2 fbc , a) bc ), 

<p ihiC = arctan 2 {a 2 ib -a] bc ). (35) 



As in the linear case, the propagator momenta p bfi are given in Eq. (27) 



D. Resonance conditions 



For the whole two-photon process, we have both momentum and energy conservation, 



as given by the four-dimensional Dirac 5 function in Eq. (26). The final electron is not 
interesting, and therefore integrated out. Left is then one constraint from the delta function. 
If this is used to fix the energy of one of the photons (we will always take photon k c to have 
fixed energy), then we are free to choose the energy u b and the direction (9 b ,ip b ) of photon 
k b , and the direction (8 C , ip c ) of photon k c . In addition, since we are interested in polarization 
resolved rates, the polarization vectors e b and e c can be chosen arbitrarily. The frequency 
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uo c can be written as a function of the direction angles 9b, 9 C , ipb, ipc as follows, 

nx ■ qi — kb ■ qi — nx ■ kf, 



nx ■ k c + qi ■ k c - k b - k c 



9 2 b E i+ ^(i+e) 



etEi + tg (i + e 



(36) 



where k c = k c /u c . In the second line of Eq. (36), we have expanded the expression for 
small io/m, 9b, 9 C and m/Ei, and we have assumed the conditions we are interested in here, 
i.e. <C |r ~ 9b ~ 9 C <C 1 (this limit corresponds to a small total exchanged laser 
photon energy as compared to the relativistic electron energy). Finally, the limiting term 
for 9b = 9 C = and Ub = is 

,max = 4nuE? 4 ni 2 u 

confirming the estimate given in Sec. | The factor (1 +^ 2 )" 1 can be interpreted simply as 
arising from the increased effective mass of the electron in the field. 

Resonances in the Dirac-Volkov propagator [TJl [38] occur if we have either p\ — m 2 = 



or pi — ml = 0. Here the two-photon amplitude (26) splits up in a product of two single 
nonlinear Compton scattering [39], HHl HI] amplitudes multiplied with a singular factor. If 
we solve for u>b, we find that the resonance conditions read 

= ~ K ' qi ~ - — ^ — , 08) 
qi -k b + sx-k b ^ + f(i+e) 

independent of n and k c (this is the usual nonlinear Compton formula |39j), and a second 
type of resonances occurs at 

res2 nx- q i -C s (q i -k c + n>c-k c ) 



a, 



-C s k b - k c + qi- kb + nx ■ k b 
sx-qi 



(39) 



qi ■ k c + s x ■ k c 



Equation (39) depends on n, so that there is one peak for each n, in principle. However, the 
dependence on s is the decisive one for typical situations. This is natural when we recall that 
s is the number of photons exchanged before the emission of the second photon. This type 
of resonance, where the electron scatters twice inside the laser pulse and emits one photon 
at each scattering event, has been referred to as "plural Compton scattering" in Ref. |42j. 



Figure [5] illustrates the formulas (38) and (39). 
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FIG. 5: (Color online.) Panel (a) shows u) r b esl [Eq. d38b] and cjf s2 [Eq. d39b], as a function of 6 b . 



Recall that we write fe& = u^(sin #6 sin sin cos ^fe, cos #;,), he. &b is the angle between k b and 
Pi. The parameters employed are ui = 2.5 eV, £ = 1, E% = 10 3 m, ipt, = ip c = 0, and 6 C = 10~ 3 . In 
panel (b), the resonance position of the first harmonic at 9 b = 1CP 3 is plotted as a function of n. 
For w£ esl , the first harmonic (the resonance at lowest possible w&) means s = 1, and since for a;£ es2 
the value of n — s tells us the order of the resonance, we have set s = n — 1 for this curve. In fact, 
for large values of n, the resonance u>l es2 with s = n — 1 shifts down to low photon energies, so that 
there will be resonances for any photon energy tJb > 0. However, these higher-order resonances will 
be suppressed by a large-order Bessel function, and effectively, one can say that the higher-order 
resonances will not contribute provided £ is not too large (~ 1). 



E. Via gauge invariance to the differential rate 



The matrix elements (26) and (31) are both invariant under the gauge transformations 

(40) 



where Ai^ are arbitrary constants (that may depend on the parameters in the problem, i.e u, 
qi etc). This symmetry can be used for a numerical check of the computer code used for the 
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evaluation, which we have performed in order to reassure ourselves regarding the consistency 
of the calculations. The gauge symmetry depends sensitively on the Bessel functions and the 
recurrence relations satisfied by them [43J, so that all signs in the formulas have to be right 
for the symmetry to hold. The gauge symmetry can also be used to simplify the expression, 
for example, by gauge transforming so that terms proportional to e& ]C • x vanish. There is 
also invariance under the transformation a — > a + Ax, A constant, but since the four- vector 



a always appears with a square, a , as ax, or in expressions like (34), this gauge symmetry 
is almost trivial and cannot be used as a meaningful validity check. 

We now discuss how to obtain the differential two-photon rates, using the example of 
linear polarization. The differential rate per unit time dW is obtained as 



dW 



T 



\S fi \ 



,Vd?q f Vd 3 k b Vd?k c 
(2tt) 3 (2tt) 3 (2tt) 3 ' 



(41) 



Here, d s kb, c = us 2 c dub, c dQb,c- The squared amplitude \Sfi\ 2 contains the Dirac 5 of argument 
zero, [o" (4) (0)] 2 = 5^( qi - q f + m - k b - k c )TV{2-n)-\ so that all factors of V and T 
in (41) cancel, as they should. We integrate over the final electron momentum and the 
photon energy u c with the delta function, and in addition we sum over the final electron 
spin (the final electron is always assumed to be unobserved), and average over the initial 
electron spin. Since in all examples we will present, the initial electron energy and laser 
intensity £ are chosen such that the quantum parameter x [ see Eq. is small, spin effects 
are marginal [35]. The final result then reads 

2 oo 4 2 2 



dW 



dwhdVLhdVLr 



X 



oo 
s=—oo 

9 9 



M s ~ n 



pi 



mi 



■x + m 



p 2 c -ml 



u n (Pi 



evaluated with qf = q% + n>c — k^ — k c and uo c is given by the first line of Eq. (36 ). Note also 
the factor QfUj c /(qf ■ k c ) arising from the delta function integration over uo c . 

In order to obtain a well-defined expression for the differential rate close to the propagator 



poles (38), (39), it is necessary to discuss some kind of regularization procedure. One 



alternative is to include an imaginary correction to the mass and energy of the laser-dressed 
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electron [HJ US] , so that Q i and m in the propagator denominator are replaced according to 

Q^Q>-i r ^f L , m ^ m -i T J^iA. (43) 

The imaginary correction is related to the total rate for nonlinear single Compton scattering 
as r(x • q) = ^Wcompton, and is given to a good approximation for small x • q/m 2 , £ = 1 
and linear laser field polarization as T(x- q) = 4 x lCT 3 x- q/m |45]. The main problem with 
this regularization scheme is that the resulting scattering amplitude is not strictly gauge 
invariant, but the noninvariance induced by the small regularizing imaginary parts of the 
energies of the virtual states is moved to higher orders. We note that very similar questions 
concerning two-photon emission amplitudes for bound states have recently been discussed in 
(S3 H3 SU HH1 [50]. As an alternative, we propose to multiply the rate with the regularizing 
factor 

oo 

$ = Yl [l - e~^( p '~ m *) ] [l - e~^( p '- m *) ] , (44) 

s=— oo 

where r is the pulse length of the laser field. The rate is now proportional to r at a 
resonance, and this way of regularizing is furthermore gauge invariant. Figure [6] shows an 



example of the differential rate (42), evaluated for a specific set of parameters corresponding 



to double Compton backscattering in a relativistically strong laser field. For cj& ^ 2 MeV, 



there is a "forest" of peaks at energies satisfying Eqs. (38) and (39). Note that according 
to Fig. [H] (b), the resonances ul es2 (s = n — 1) should actually lead to resonance peaks also 
at very low photon energies, but these are suppressed by large-order Bessel functions and 
thus not visible. The bright curves in Fig. [6] (c) correspond to the maxima in the differential 
rate induced by single-Compton scattering, but the rate is nonvanishing in other areas of 
the ^-co^-plane due to the two-photon emission. 

The object of this paper is however not to study the behavior of the process close to the 
peaks, but rather to single out a kinematic region where unambiguous conclusions can be 
drawn independent of the method of regularization. The kinematic region best suited for 
such investigations seems to be for photon energies Ub and angles smaller than some thresh- 
old such that the contribution from the cascade peaks are negligibly small. Mathematically, 
the suppression arises due to a large-order generalized Bessel functions (or, alternatively, 
ordinary Bessel functions in the case of circular laser polarization), which beyond some cut- 
off index decays exponentially with increasing n [37| |4~3] . In all subsequent examples in the 
remaining sections of this paper, we will therefore restrict the photon energy Uj, and the 

16 



50 100 
u. b [MeV] 



150 



(b) 



10 



^10 M - ~ 



d io 

:: 

J 10" 



Meth. 1 




Meth. 2 


1 * 

ii H i 


j jJJ 





2.3 2.35 2.4 2.45 

Ul h [MeV] 




0.5 1 1.5 2 2.5 3 
"° L0 b [MeV] 



3.5 




FIG. 6: (Color online.) We illustrate the resonances of the fully differential rate (42) as a function 
of cob. The parameters used are £ = 1, Ei = 10 3 m, = e c = ej, 6 C = 0.002, ipf, = ip c = 0. In 
panels (a) and (b), a fixed value of 9b = 0.001 is used, corresponding to a cut along the dashed line 
in (c). In panel (c), the color coding indicates the value of the decadic logarithm log 10 dw b dn b dn ' 
where the argument of the logarithm is measured in units of s~ 1 sr~ 2 MeV _1 . The regularization 
method employed is given by a finite laser pulse duration as in Eq. d44l, with r = 10 4 /w. In the 



inset of panel (b), both methods (43) [Meth. 1] and ( |44[ ) [Meth. 2] are shown for comparison. For 
ujb < 2 MeV, the application of the two methods yields numerically indistinguishable curves. The 
first Compton harmonic [the lowest bright curve in panel (c)] is broken at 9b ~ 1-2 x 10 -3 , which 
can be understood as the point where a ■ €b = in the rest frame of the electron [in this frame, 
and in the gauge eb = (0, e^), the Thomson cross section is proportional to \a ■ eb\ 2 ]. 
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polar angle 0& to the region u b < 1 MeV and 6*b < 0.002. Here, the result is independent of 
the method of regularization since we are sufficiently far away from the cascade peaks. With 
increasing £, the "safe" region shrinks, as the first Compton peak appears at lower energy 



LOb, see Eq. (38). Already at £ = 2, there are cascade contributions at tu b < 1 MeV, why we 



limit the laser intensity to £ < 1 in the following. 



III. COMPARISON TO PERTURB ATIVE DOUBLE COMPTON SCATTERING 



In the limit £ — > 0, the amplitudes (26), (31) reduce to the one found in [5], where only 



one photon is absorbed from the laser. A discussion of this process can be found in standard 
textbooks [51], and was recently reexamined in [3]. The above mentioned references as 
well as other previous works [TT| l52| 153] were devoted to the study of the cross section for 
unpolarized initial and final photons (with a few exceptions, see j5U|55]). However, as noted 
in [18], the discussion of photon polarization correlation necessitates an expression for the 
amplitude for arbitrary polarization of the final photons. 

The amplitude S PDCS for perturbative double Compton scattering (PDCS) is given by 
the sum of the Feynman diagrams shown in Fig. [TJ and reads [51] 



s PDCS = ™t 3(27r)4 ,J (Pf + h + k c -x- P . 



^%E i EfUJUJ b uj c V t 

x>y (45) 



.1=1 



where 



. r _/ \ a p f + k c + m . pi + k + m 

M -i \ ~ Pf + k + m A Pi -k b + m A . . 

N2 = U{Pf)£c ( Pf + k c )i- m i\ Pl -k b r- m i ebU{P ^ (46) 
, r _/ ^ pf-x + m , pi-k b + m 

N * = u ^ ) \ P} -xy- m ^{ Pl -k b y- m * ebU{Pi) ' 

A 4 = Ni(b <->• c), N$ = N 2 (b <-» c) and A 6 = N 3 (b <-> c) in a self-explanatory notation. Here, 

e represents the polarization vector of the laser field. The rate is then 

7 I cPgC5|2 1 , r Vd 3 h Vd 3 k c Vd 3 Pf 
aW ~ J \ b I T V ( 2vr )3 (2vr)3 (2vr)3 



4„4 



4 me 



8EiE fU b u c (2ir) 



i=i 



<KP/ + ^6 + — Pi — x)d 3 k b d 3 k c d 3 P f, 
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where J = is the photon flux. Integrating over d 3 pf and ckd c , and summing and 

averaging over the electron spin, we obtain 

(48) 



du b dtt h dtt c ~ * 16(27r) 5 Ei)o f • k c ^ 

v ' J r;,r/=l 



The dependence on £ is thus given by the prefactor £ 2 , and we may thus refer to this process 
as the "linear" process because the rate is proportional to the laser intensity. 

In Figs. [7] and [8j we show a comparison of the predictions of the nonperturbative formulas 



[Eq. (26) and Eq. (31)] and the perturbative formula [(45)] for the differential rate, for both 
circular and linear polarization of the laser field. Note that to compare with the circular 
polarization, one should put e = (Q,l,i,0)/y/2 in Eq. ((461. The fi gures show that the 



photons may be produced in any polarization state, although parallel polarization of the 
two emitted photons e& • e c ~ 1 is the dominant channel for linear laser polarization. Due 
to the rotational symmetry, circular laser polarization results in similar differential rates 
for both parallel and perpendicular polarization of the emitted high-energy photons [see 
Figs.[7](d)-(f) and (j)-(l)]. For small laser intensity, the plane defined by the polarization and 
propagation axes of the laser characterizes the emission pattern of the accelerated charge, 
so that alignment of the two emitted photons as shown in the third row of Fig. [7] can 
be expected. Our numerical results show that this intuitive picture is still valid in the 
relativistic, nonperturbative laser interaction regime, although some details change: e.g., 
the emission of photons with antiparallel polarization vectors is favored over the parallel 
case, as is evident from the first panel in the upper row of Fig. [7j We here recall that the 
direction of the polarization basis vectors e l b ,2 c depends on the angles ipb tC , see Fig. 3j Figure ^ 
confirms this picture, here the difference of the differential rate between the parallel case 
(e& • e c ~ 1) and the perpendicular case (e& • e c « 0) amounts to several orders of magnitude. 

In order to demonstrate the contribution from the different photon orders n, we refer to 
Fig. [9j Here the dependence on the photon order n is shown, if we define 

dw =y dwn ( 49) 

diObdVLbdVlc 2 -— ' dujbdilbd^lc 

n=l 

As can be seen from Fig. |9l typically up to 20 photons contribute to the differential rate. 



Equation (36) for n < 20 then yields uo c < 70 MeV, which implies that even though the 



"first" photon has modest energy uj, — 1 MeV, the energy of the "second" photon is much 
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larger. The difference between the smooth curve in the circular case and the sawtooth shape 
of the linear curve can be traced back to the behavior of the generalized Bessel function and 



the usual Bessel function, constituting the amplitudes (26) and (31). For example, for the 



parameters shown in Fig. |9[b), one can show that the dominant contribution to the matrix 
element for linear polarization is roughly proportional to the generalized Bessel function 
A±(n, 0, (3), which vanishes for even n [13]. However, if the polarization vectors are summed 
over, then the case of even n contributes, and the curve smoothens out. Similar selection 
rules for the emitted harmonics occur also for the nonlinear single Compton scattering 
process [35J . On the contrary, the circular polarization curve is smooth due to the rotational 
symmetry. 

To conclude this section, we investigate if the integrated rate differ in the perturbative 



and nonperturbative case. In Fig. [TOj we show, as a function of £, the quantity 

27T 27T ^b,max 



W int = ^2 J d^b J dip c J sin 6 b d6 b 

^hi^n n n n 

(50) 







sin 9 c d9 c I dcu b 



duJbdVtbdVL c 

<^6,min 

with 6» bjmax = 1.5 x 10" 3 , 6> c>max = 2.5 x 10" 3 , u bMn = 10" 3 MeV, and u; bjmax = 1 MeV. 
This restriction is identical to the one in [31]. By restricting the final phase space, one 
can ensure that contributions from the single Compton scattering cascade are negligible, as 



discussed in Sec. |II E| Figure [10] reveals that the integrated rate is slightly larger than one 
would expect from the perturbative formula, and also that circular and linear polarization 
of the laser gives almost identical results, despite their different angular characteristics (see 
Figs. [7j [8]). Another remark is that for the integrated rate, the perturbative formula gives 
identical results regardless of laser polarization because interference terms in the expanded 
perturbative rate vanish after the integration. In the interval considered (£ < 1), the relative 
difference of the integrated nonperturbative and the integrated perturbative rate can be 
approximately fitted to a power law as W^ npevt — W^ vt oc with rj ~ 2.7 for linear and 
i] m 3 for circular laser polarization, respectively. 
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FIG. 7: (Color online.) Contour plot of the differential rate for different final photon polarizations. 
The color coding indicates the value of the differential rate g^~gfy~gQ i on linear scale, in units of 
s _1 sr _2 MeV~ 1 . The upper two rows [(a)— (f)] display the result of the nonperturbative formula, and 
the lower two rows [(g)— (1)] show the perturbative results. The parameters used are £ = 1, u = 2.5 
eV, Ei = 10 3 m, uj^ = 1 MeV, and 9b = 26 c = 10 -3 . In this case, uj c is fixed by the scattering 
geometry but still depends on the number n of exchanged photons; we present the differential 
rate summed over all n and thus suppose that the energy u c is unobserved. For the polarizations, 
we have linear laser polarization in (a)-(c) and (g)-(i) (first and third row), and circular laser 
polarization in (d)-(f) and (j)-(l) (second and fourth row). In (a), (d), (g), and (j) (first column) 
we have et, = e^, e c = e\, where we recall that the polarization vectors are defined in Eq. ( |24[ ). In 
(b), (e), (h), and (k) (middle column) we have e& = e\, e c = e 2 , and panels (c), (f), (i) and (1) (right 
column) show the differential rate summed over all possible polarizations of the final photons. 
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FIG. 8: (Color online.) Contour plot of the differential rate as a function of the polar angles 9^ 
and 6 C . The color coding indicates the value of the decadic logarithm log 10 d^ b dQ b dn > w hh the 
differential rate given in units of s -1 sr — 2 MeV~' 1 '. As in Fig. [7J the upper two rows [(a)-(f)] show 
nonperturbative, and the lower two rows [(g)— (1)] perturbative results, respectively. The parameters 
used are £ = 1, u = 2.5 eV, E{ = 10 3 m, ut, = 1 MeV, and tpb = 0, ip c = it. For the polarizations, 
we have linear laser polarization in (a)-(c) and (g)-(i) (first and third row), and circular laser 
polarization in (d)-(e) and (j)-(l) (second and fourth row). In (a), (d), (g), and (j) (left column) 
we have e;, = e^, e c = ej, In (b), (e), (h), and (k) (middle column) we have e;, = e^, e c = e 2 , and 
panels (c), (f), (i) and (1) (right column) show the differential rate summed over the final photon 
polarizations. Note that from Fig. hi we have e\ ■ e\ ~ 1 and • e 2 ~ here. 



IV. ANGULAR CORRELATION AND ENTANGLEMENT 

We now turn our attention to the important questions regarding the quantum mechanical 
correlation, i.e. entanglement, of the two final photons. The theory we apply in this section 
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FIG. 9: (Color online.) Demonstration of the multiphoton character of the pair creation rate, for 
linear and circular polarization of the laser field. In (a), tpb = 0, ip c = and the final polarization 
is fixed to €b = e c = ej, while in (b) we have ipb = n/2, ip c = 3n/2, and = e^, e c = e 2 c . 
Otherwise the parameters are the same as in Fig. [7} 



have been previously used extensively to characterize the final-state correlation in bound 
states transitions [HSl E3, EEJ EH] • The idea is to use the information contained in the matrix 



elements (26), (31), to obtain an expression for the density matrix pf of the polarizations 
of the final system "electron+two photons". Given an expression for pf, it is then straight- 
forward to calculate the concurrence [19J, which is a measurement of how much the two 
photons are entangled. The starting point is the initial density matrix [60] 

2 

^ = ^1^,0,0)^,0,01, (51) 
n=i 
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FIG. 10: (Color online.) The integrated rate (50 ), denned in Eq. (50 ) for lj = 2.5 eV and Ei = 10 m. 



The figure compares the nonperturbative formula for linear laser polarization [Eq. (42 )], for circular 



laser polarization [Eq. (31)], and the perturbative expression [Eq. ([48])] . 



where is the spin of the initial electron and the zeros denote the absence of photons (other 
than laser photons of course) in the initial state. The initial electron is thus assumed to be 
unpolarized. Note also that all dependencies on energies and angles etc of the state vectors 
are not written out. Next, due to the interaction R, the density matrix pi evolves into the 
final state density matrix pf, 

2 

Pf = R Pl rf = R K 0, 0)(ru 0, 0\Rl (52) 

The matrix elements of pf are thus given by 

(77, X b ,X c \p f \r' f , X' b , X' c ) 
2 

where Xb, c , X' bc G {1,2} denotes the polarization components of the emitted photons in either 
Cartesian or circular basis. If the final electron is unobserved, we should trace out rf. 

(x b , \c\pAKiK) 
2 

= (rf,X b ,Xc\R\ri,0,0)(r u 0,0\R^\r f ,X' b ,X' c ). (54) 
n,rf=i 
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If we now identify 

(77, A 6 , \ c \R\r u 0,0) = y/NSfi(ri,rf,Xb,\c), 
where N is a normalization constant, and we use the explicit basis 



|H) 



|12) 



|2 1> = 



|2 2) 






/o\ 

1 


/o\ 


1 

/o\ 




w 



(55) 



(56) 



for the polarization state of the final photons, then the expression for the final density 4x4 
matrix reads 

/ c c* go* go* c c* \ 

OnD 11 OiiD 12 <JHD 21 »J11>J22 

go* go* go* go* 

Oi2<J n Oi2>J 1 2 01202! Oi20 22 

C C* c 0* 00* 00* 

021>J 11 021>J 1 2 >J21'- , 21 >- , 21*J22 

c c* c c* go* c c* 

\ <J22>Jll >J22>Jl2 >J22*->21 >- , 22>->22 / 



p/= E N 

n,rf=l 



(57) 



where 



S\ b \ c = S fi (ri,rf,\ b ,\ c ). 
The normalization constant iV can be found by requiring 

1 = Tr P/ = N (l^iil 2 + ISd 2 + |^2i| 2 + |^ 22 | 2 ) 



(58) 



(59) 



Ti,r f 
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According to [19], the concurrence C{pf) is now given by 

C(p f ) = max(0, a/Ci - a/G - VC3 - y/U), 
where the Q's are the eigenvalues, in descending order, of the matrix 

Q = pf {<y 2 ® (x 2 ) p) (a 2 <g> a 2 ) , 



(60) 



(61) 



where a 2 is the second Pauli spin matrix. The eigenvalues Q are real and positive. The 



concurrence as defined in Eq. (60) is gauge invariant, and does not depend on the basis used 



for the polarization vectors of the photons, i.e. either the Cartesian basis, Eq. (24), or the 
helicity basis, Eq. (25), can be used. An explicit expression for (a 2 ® a 2 ) as a 4 x 4 matrix 
is given by 

/ -1 \ 
10 
10 
\-1 / 

The matrix (a 2 <8> <J 2 ) is a kind of spin-flip operator for qubits, we have 



( 



a 2 ® a 2 ) 







fo-i\ 




M 




\< 01 




I i j 



(62) 



( 



a ® a 



') |1 1) = -|2 2), [a 2 ® a 2 ) |2 2) = -|1 1), 



[a 2 ® a 2 ) |1 2) = |2 1), [a 2 ® a 2 ) |2 1) = |1 2). 
This means that a maximally entangled pure state is an eigenstate of (cr 2 <g> a 2 

(a 2 ® a 2 ) (|1 2) - |2 1)) = -|1 2) + |2 1), 

and consequently has unity concurrence. 



(63) 



(64) 



We now provide some examples and compare the concurrence (60) for different laser 



polarizations and furthermore show that the nonperturbative treatment is indispensable to 



correctly predict the degree of correlation. Figure 11 shows the concurrence as a function 
of the azimuth angles ipb, ip c . This figure should be compared to Fig. [7j In Fig. [T2j we 
show instead the dependence on the polar angles 9^ c , which should be compared with the 
corresponding Fig. [8] for the differential rate. 

We remark that to be able to measure the concurrence, it is desirable to find angular 
regions where high concurrence and high differential rate overlap. This seems to be possible, 
at least in some cases: e.g., one may compare Fig.[l27c) with Fig.[8tl). Moreover, the general 
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trend is that a strong laser field diminishes the concurrence. Therefore, if high entanglement 
is sought, it is advisable to employ a perturbative laser beam (£ < 1), although the non- 
perturbative dependence of the concurrence as a function of £ would be highly interesting 
to measure. Similar conclusions as those above follow from the previous investigation |34j. 
A final remark is that linear and circular laser polarization are seen to lead to similar peak 
values of the concurrence. 




FIG. 11: (Color online.) Contour plot of concurrence as a function of the azimuth angles ipb and ip c . 



The color coding indicates the value of the concurrence C(pf) [see Eq. (60)]. Panel (a) shows linear 
laser polarization, and (b) circular laser polarization, both calculated with the nonperturbative 
expressions for the amplitude. The results of the perturbative formula are displayed in panel (c) 
[linear laser polarization] and (d) [circular laser polarization]. The values of Ei, to, £ and 9^ c are 
the same as in Fig. [7] 



V. CONCLUSIONS 

Our treatment of double Compton scattering in intense laser fields in based on the canon- 
ical formalism of Furry-picture quantum electrodynamics (QED), where a strong external 
field (in this case, the oscillatory laser field) is incorporated into the fermion propagators. 
The oscillatory nature of the laser field necessitates the expansion of all initial and final 
fermion states into plane waves, thereby giving rise to generalized (linear polarization) and 
ordinary (circular polarization) Bessel functions. The formalism, as outlined in Sec. [TlJ leads 
to a consistent formulation of the nonperturbative double Compton scattering for an arbi- 
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FIG. 12: (Color online.) Contour plot of concurrence as a function of the polar angles 9b and 9 C . 



The color coding indicates the value of the concurrence C(pf) [see Eq. (60)]. Panel (a): linear laser 
polarization, nonperturbative, (b): circular laser polarization, nonperturbative, (c): linear laser 
polarization, perturbative, (d): circular laser polarization, perturbative. Laser parameters as well 
as the initial electron energy Ei and the photon azimuth angles ipt>,c are the same as in Fig. |8j 



trary intensity of the laser field, and with full account of all relativistic and spin-dependent 
effects on the electron lines. In particular, a suitable generalization of the formalism outlined 
here would apply to three-photon events, which can be described by a third-order amplitude 
in QED. 

In addition to a consistent formulation of the polarization resolved production rates, 
differential in the photon emission angles and energy, for two-photon transitions of Dirac- 
Volkov states in intense laser fields, we numerically show that only a fully relativistic formal- 
ism, nonperturbative in the laser field strength, can possibly yield experimentally verifiable, 
consistent predictions. This is not surprising because the differential rates depend crucially 
on details of the emission process, which in turn is highly dependent on the properties of the 
propagators near the resonances. Indeed, it is possible to identify those angular and photon 
energy regions where only the two-photon amplitude, not the single-Compton resonances, 
give appreciable contributions to the photon emission, and it is thus possible to observe 
entangled high-energy photons in coincidence without having to worry about background 
from resonant cascade emission by single-photon transitions. 



The necessity of the nonperturbative formalism is demonstrated in Sec. |III[ The per- 
turbative (in the laser field interaction) double Compton scattering cannot give reliable 
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predictions if the nonlinear intensity parameter £ approaches unity. We stress that a value 
£ ~ 1, corresponding to laser intensities of the order of 10 18 W/cm 2 for optical lasers, is 
routinely available today in many laboratories worldwide. In Figs. [7] and |8j we show that 
depending on polarization and the observation solid angle, even order-of-magnitude differ- 
ences can exist between the rates evaluated with the nonperturbative and the perturbative 
formulas. However, if angles and photon energies are integrated over, the results are simi- 



lar, as shown in Fig. [TOj although the difference grows nonlinearly with £, illustrating the 
importance of higher orders. 

The polarization entanglement is interesting but needs to be quantified. Therefore, we 



discuss, in Sec. IV, the concurrence as a gauge-independent measure of the photon en- 



tanglement. Our results (see Figs. 11 and 12) indicate that close to maximally entangled 



(unity concurrence) photon pairs may be produced, but only in certain angular regions. 
Furthermore, the degree of entanglement changes strongly with the laser field intensity. 
An experimental verification of the entanglement would yield a test for this fundamental 
quantum phenomenon in a high-energy domain where it is otherwise difficult to generate 
entangled quanta. 

Finally, we remark that the two-photon emission is not a "rare" or "unusual" physical 
process but a simple generalization of the basic physical phenomenon of radiation emission 
by moving charges, and that, therefore, we can assume that experimental access in the near 
future is entirely realistic. 
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